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We use maximally localized Wannier functions to construct tight-binding (TB) parameterizations 
for the Cg bands of LaMnOs based on first principles electronic structure calculations. We compare 
two different ways to represent the relevant bands around the Fermi level: i) a d-p model that includes 
atomic- like orbitals corresponding to both Mn{d) and 0{p) states in the TB basis, and ii) an effective 
Cg model that includes only two e^-like Wannier functions per Mn site. We first establish the effect 
of the Jahn- Teller distortion within the d-p model, and then compare the TB representations for 
both models obtained from GGA+[/ calculations with different values of the Hubbard parameter U . 
We find that in the case of the d-p model the TB parameters are rather independent on the specific 
value of [/, if compared with the mean- field approximation of an appropriate multi-band Hubbard 
Hamiltonian. In contrast, the U dependence of the TB parameters for the effective Cg model cannot 
easily be related to a corresponding mean-field Hubbard model, and therefore these parameters 
depend critically on the specific value of [/, and more generally on the specific exchange-correlation 
functional, used in the electronic structure calculation. 



I. INTRODUCTION 

The construction of realistic low-energy Hamiltonians 
based on first principles electronic structure calculations 
is an important tool for the investigation of correlated 
electron systems, such as e.g. cuprates, manganites, or 
other complex transition metal (TM) oxides. While much 
of the general physics determining the diverse properties 
of these materials can be understood in terms of simpli- 
fied models the corresponding model parameters are 
generally unknown and have to be adjusted to fit ex- 
perimental data. Alternatively, electronic structure cal- 
culations based on density functional theory (DFT)I^ 
can be used to determine these parameters,!^ which in 
turn allows for a quantitative evaluation of the underly- 
ing model assumptions and the construction of realistic 
model Hamiltonians with all parameters determined ah 

Model Hamiltonians for correlated electron systems are 
typically formulated within a tight-binding (TB) picture 
that involves only a small number of electronic states 
localized on certain atoms (e.g. the "(i-states" of the 
TM atoms within a TM oxide) .EH^ Ideally, these states 
give rise to an isolated group of bands around the Fermi 
energy, which then determines the low energy behavior 
of the system. A corresponding TB representation can 
in principle be obtained by constructing Wannier func- 
tions from the Kohn-Sham Bloch states calculated within 

In the case of model Hamiltonians that contain an ex- 
plicit electron-electron interaction, typically in the form 
of a local Hubbard term with interaction parameter /7, 
the DFT band-structure can either be viewed as mean- 
field approximation to this interacting model, or as repre- 
sentative for the "non-interacting" case, i.e. correspond- 
ing ioU = in the model Hamiltonian. While it is often 
convenient to consider the electronic structure calculated 
within either the local density approximation (LDA)P^ 



or the generalized gradient approximation (GGA^P^ as 
essentially non-interacting, this is probably not a good 
assumption if the electronic structure is calculated using 
more advanced exchange correlation functionals such as 
(LDA/GGA)+/7 or hybrid functionals .E^HSI 

In this work we use maximally localized Wannier func- 
tions (MLWFs)^^ to obtain TB parameterizations for 
the important case of LaMnOa, the parent material of the 
colossal magneto-resistive manganites, and a prototype 
system for correlated electron physics. f^^^ We calculate 
the electronic structure of LaMnOa using the GGA+/7 
method and different values for the Hubbard parameter 
U . We then compare two different ways to represent the 
relevant bands around the Fermi level: i) a d-p model 
that includes atomic-like orbitals corresponding to both 
Mn((i) and 0{p) states in the TB basis, and ii) an ef- 
fective Cg model that includes only two e^-like Wannier 
functions per Mn site. In particular, we analyze the U 
dependence of the two different TB parameterizations 
and relate this to commonly used model Hamiltonians 
for manganites. 

We find that for the d-p parameterization, the effect of 
U is mostly local, leading to U dependent shifts of the 
on-site energies and an increase in the Jahn- Teller (JT) 
splitting, whereas the corresponding hopping amplitudes 
are only weakly affected by the value of U. Thus, the 
U dependence of the d-p parameters closely resembles 
the U dependence of a corresponding multi-band Hub- 
bard Hamiltonian in mean-field approximation. On the 
other hand, the change of on-site energies and JT split- 
ting for the effective Cg MLWFs, calculated for different 
values of [/, are distinctly different from a correspond- 
ing mean-field model Hamiltonian. In addition, there is 
also a strong U dependence of the effective Cg hopping 
amplitudes, which is due to electronic degrees of freedom 
that are excluded from the TB basis. This indicates that 
while the TB parameterization for the d-p model is fairly 
robust with respect to a variation of /7, for the case of 
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the effective Cg basis an appropriate choice of U in the 
GGA-\-U calculation is crucial. Our results demonstrate 
the simple fact that the transferability of a specific TB 
parameterization is generally increased if more electronic 
degrees of freedom are included in the model description. 

This paper is organized as follows. In the next section 
we first summarize the most important aspects of typi- 
cal model Hamiltonians for manganites, emphasizing in 
particular the role of the J T dis tortion (Sec. II A), before 
introducing MLWFs (Sec. II B) and describing the tech- 



nical details of our calculations (Sec. II C). In Sees. Ill A 



and IIIB we first present the main features of the d-p 
model based on the MLWFs calculated for U = eV, 
and then clarify the effect of the JT distortion on the 
corresponding TB parameterization. An analogous anal- 
ysis for the effective Cg model parameterization has re- 
cently been presented in Ref. HI The effect of varying the 
Hubbard U in the GGA+[/ calculation on the MLWF 
parameters of the two different models is presented in 



Sec. IIIC for the d-p model and in Sec. HID for the ef- 
fective Co model. Finally, Sec. llV 



summarizes our mam 



conclusions, while further details about the d-p TB pa- 
rameterization can be found in Appendix [A) 



II. METHOD AND THEORETICAL 
BACKGROUND 

A. Jahn- Teller distortion in LaMnOs 

LaMnOs crystallizes in an orthorhombically distorted 
perovskite structure with Pbnm space group symmetry.^ 
The distortion relative to the cubic perovskite structure 
can be decomposed into three components P a staggered 
JT distortion of the MnOe octahedra, alternating tilts 
and rotations of these octahedra around the cubic axes 
(GdFeOs distortion), and an orthorhombic strain of the 
unit cell. The electronic structure of LaMnOs around the 
Fermi energy is dominated by Mn 3d states, which are 
split by the octahedral c rystal field into lower-lying t2g 
and higher-lying Cg states. The formal occupation 
of the Mn^+ cation leads to a high spin configuration 
with filled majority spin t2g states, half-filled majority 
spin eg states, and empty minority spin states. The low 
energy behavior of LaMnOs is therefore governed by the 
partially filled majority spin bands with predominant Cg 
character. 

Motivated by this, the complex phenomenology ob- 
served in manganites is often modeled within an effec- 
tive TB model that involves only two Cg orbitals per 
Mn site. Electrons can then hop between the two Cg 
levels on neighboring sites, and can interact with each 
other through a Hubbard-type electron-electron interac- 
tion, with the t2g "core-spins" through a Hund's rule 
interaction, and with the local JT distortion through a 
crystal field splitting (see e.g. Ref. 3). It is usually under- 
stood that the corresponding "e^ orbitals" are spatially 
extended Wannier orbitals that result from hybridization 








\ 















FIG. 1: (Color online) Left: ideal cubic perovskite structure 
viewed along the z axis. Right: staggered Q'^-type JT dis- 
tortion within the x-y plane. Oxygen anions are shown as 
small (red) spheres, the La cation as large (green) sphere. 
Mn cations (not shown) are situated in the middle of each 
oxygen octahedron. 



between atomic-like Mn(e^) orbitals and the p orbitals of 
the surrounding oxygen ligands. 

The effect of the JT distortion on this Cg manifold is 
typically expressed via a local crystal-field splitting of the 
form: 
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Here, A is the JT coupling strength. Tab are the usual 
Pauli matrices, ci^Yia is the annihilation operator corre- 
sponding to orbital b with spin a at site R, and the JT 
distortion of the oxygen octahedron surrounding site R 
is described by specific modes QJ^ {i = x,z)W 

In the following we only consider a staggered JT dis- 
tortion of the form = — Q^/, where R and R' corre- 
spond to nearest neighbor sites within the x-y plane (see 
Fig.[l]). Thereby: 
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with d^ and being the Mn-0 distances along the x 
and y directions, respectively, corresponding to the oxy- 
gen octahedron surrounding the Mn at site R. We note 
that it has been shown in Ref. 8 that this component of 
the JT distortion has the most pronounced effect on the 
electronic structure of LaMnOs. Furthermore, in Ref. [9] 
we showed that the effect of the various structural distor- 
tions on the calculated model parameters can be analyzed 
separately, since they are to a great extent independent 
of each other. 



B. Maximally localized Wannier functions 

A set of N localized Wannier functions \wnT) corre- 
sponding to a group of N bands that are described by 
delocalized Bloch states IV^mk), is defined by the follow- 
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ing transformation: 
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(3) 



Thereby, T is the lattice vector of the unit ceh associ- 
ated with the Wannier function, m is a band index, k 
is the wave-vector of the Bloch function, and the inte- 
gration is performed over the first Brillouin zone (BZ) of 
the lattice. Different choices for the unitary matrix U^^^ 
lead to different Wannier functions, which are thus not 
uniquely defined by Eq. ([3|. A unique set of maximally 
localized Wannier functions (MLWFs) can be generated 
by minimizing the total quadratic spread of the Wannier 
orbitalsP 

Once the transformation matrices U^^^ are deter- 
mined, a TB representation of the Hamiltonian in the 
MLWF basis is obtained: 
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Here, e^k is the eigenvalue corresponding to Bloch func- 
tion I'^/^^k)- For cases where the bands of interest do not 
form an isolated set of bands but are entangled with other 
bands, a two step procedure for obtaining the unitary 
transformation matrices (which in this case are typically 
rectangular) is employed. 

We note that T and AT in Eqs. indicate lat- 

tice translations, whereas for crystal structures with more 
than one atom per unit cell n and m represent a combined 
orbital and site index, specifying the various orbitals at 
all sites within the primitive unit cell. 



C. Computational details 

All results presented in this work are obtained 
from spin-polarized first principles DFT calculations us- 
ing the Quantum-ESPRESSO program package, the 
GGA exchange-correlation functional of Perdew, Burke, 
and Ernzerhof,^^ and Vanderbilt ultrasoft pseudopoten- 
tials.^La (55, 5p) and Mn (35, 3p) semicore states are in- 
cluded in the valence. The Hubbard -\-U correction is 
applied using the simplified approach according to Du- 
darev et which corresponds to the case J = in 

the more elaborate expression by Lichtenstein et al.^ 
Projections on orthogonalized atomic Mn{d) orbitals are 
used to evaluate the U dependent contributions to po- 
tential and energy. 

To analyze the effect of the JT distortion on the elec- 
tronic structure of LaMnOa, we perform GGA calcula- 
tions with different degrees of distortion. Starting from 



the ideal cubic perovskite structure, we gradually in- 
crease the amplitude of the JT distortion \Qji\ from 
to Qq = 0.151 A. The latter value corresponds to the 
amount of distortion found in the experimentally ob- 
served crystal structure of LaMnOa.^^ For the undis- 
torted case we use a cubic perovskite structure with lat- 
tice constant ao = 3.9345 A, which results in the same 
volume V = 60.91 A^ per formula unit as in the experi- 
mentally observed Pbnm structure.^ 

After obtaining the DFT Bloch bands, we construct 
MLWFs using the wannier 90 program integrated into 
the Quantum-ESPRESSO package. Starting from an 
initial projection of atomic d and p basis functions cen- 
tered on the Mn and O sites onto the Bloch bands within 
an appropriately chosen energy window, we obtain sets 
of either 14 atomic-like {d-p model) or 2 more extended 
e^-like (effective Cg model) MLWFs per spin channel and 
unit cell. In both cases we will in the following refer to the 
MLWF Hamiltonian matrix elements between the 
two different e^-like MLWFs located on the same Mn site 
as on-site off-diagonal element to any matrix element 
between the same orbital at a particular site as on-site 
energy £, and to matrix elements connecting two differ- 
ent sites as hopping amplitudes t. To assess the effect of 
the correction on the two different MLWF parame- 

terizations we then construct the corresponding MLWFs 
for the fully JT distorted structure from GGA+/7 calcu- 
lations with different values for the Hubbard U. 

Convergence of the DFT total energy and total mag- 
netization has been tested for the ideal cubic perovskite 
structure and ferromagnetic (FM) order. We find the to- 
tal energy converged to an accuracy better than 1 mRy 
and the total magnetization converged to an accuracy 
of 0.05 /iB for a plane- wave energy cut-off of 35 Ry and 
a F-centered 10 x 10 x 10 k-point grid using a Gaussian 
broadening of 0.01 Ry. These values for plane- wave cut- 
off and Gaussian broadening are used throughout this 
work, whereas an appropriately reduced k-point grid of 
7x7x10 is used for the JT distorted structure with 
doubled unit cell in the x-y plane. The MLWFs are 
considered to be converged if the fractional change of 
the quadratic spread (both gauge-invariant and non- 
gauge- invariant part) between two successive iterations 
is smaller than 10~^^. 



III. RESULTS AND DISCUSSION 

A. d-p TB parameterization for cubic LaMnOa 

In this section we establish the general features of the 
extended d-p TB description of LaMnOs within the ideal 
cubic perovskite structure, before we analyze the effect 
of the JT distortion in the next section. We are consider- 
ing a FM arrangement of magnetic moments, but we have 
verified that the corresponding results for A-type antifer- 
romagnetic (A-AFM) order do not exhibit any significant 
differences. 
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FIG. 2: (Color online) Projected DOS and band structure 
along high symmetry lines within the BZ for cubic FM 
LaMnOs. In the projected DOS plots, the filled (red) ar- 
eas and solid (green) lines correspond to Mn{eg) and Mn{t2g) 
states, respectively, while dotted (blue) and dashed (orange) 
lines correspond to the O(pcr) and 0(p7r) states, respectively. 
In the band structure plots, the dispersion calculated from 
the MLWFs is represented by thick (red) lines. The Fermi 
level is indicated by the horizontal dashed lines. 



Fig. [2] shows the projected densities of states (DOS) 
and band dispersion for both majority and minority spin 
channels. In agreement with previous calculations it can 
be seen that the majority spin DOS around the Fermi 
energy exhibit predominant Mn{eg) orbital character, 
with M n(t2Q) an d 0{p) bands located at slightly lower 
energiesj ^ l ^ l ^^ l ^^ l For minority spin, states with predom- 
inant Mn(e^) and Mn{t2g) character are located above 
the Fermi level. Strong hybridization between 0{p) and 
Mn{d) orbitals is apparent from the various peaks in the 
DOS around 8 eV. 

We construct 14 MLWFs per spin channel from the 
Kohn-Sham states located within an energy window of 
3-17 eV and 3-20 eV for majority and minority spin, 
respectively. This corresponds to a TB representation 
of LaMnOs containing 5 Mn{d) orbitals and 9 0{p) or- 
bitals per unit cell (three p orbitals corresponding to each 
of the three oxygen atoms within the simple cubic per- 
ovskite unit cell). The resulting MLWFs are depicted 
in Fig. [3] Three of the 0{p) orbitals (one on each O 
atom) are pointing towards the Mn atom and hybridize 
with the Mn(e^) orbitals. We call them 0{pa) orbitals. 
The remaining 6 0{p) orbitals (two on each O atom) are 
oriented perpendicular to the Mn-0 bond and hybridize 
with the Mn{t2g) orbitals. These will be called 0(^7^) 
orbitals in the following. 

The MLWF bands (shown as thick (red) lines in Fig. [2| 
are identical to the corresponding DFT bands for ma- 
jority spin, whereas for minority spin the MLWF and 
DFT bands above the Fermi energy exhibit certain dif- 
ferences which are due to the strong entanglement with 
other bands in that energy region. 

From the real-space Hamiltonian matrix elements in 
the MLWF basis, Eq. ([5|, we find that the most domi- 





FIG. 3: (Color online) Real space representation of the ML- 
WFs corresponding to the d-p parameterization (FM, major- 
ity spin). Depicted are Mn(32:^— r^), Mn(x^— 2/^), and Mn(x2/) 
MLWFs (top row from left to right) as well as two examples 
of O(pcr) and 0(p7r) orbitals (bottom row). The isosurface 
shown corresponds to a value of ^Xj^/V where V is the unit 
cell volume. Picture generated using XCRYSDEN.^^- 



nant hopping corresponds to the shortest Mn-0 and 0-0 
bonds (see Fig. 10), but that a variety of other hopping 
amplitudes are also non- negligible. In the following we 
will focus on the nearest neighbor Mn-0 hopping am- 
plitudes {teg, pa and tt2g,p7T) ^ud analyze how they are 
affected by the JT distortion and the inclusion of a Hub- 
bard U in the DFT calculations. A more detailed discus- 
sion of the d-p TB parameterization for cubic LaMnOs 
can be found in Appendix [Aj 



B. Effect of JT distortion on the d-p TB 
parameter izat ion 

We now analyze the effect of the JT distortion on the 
on-site energies and nearest neighbor hopping parame- 
ters of the d-p model. As in the previous section, we 
will discuss only results for FM order as we do not find 
significant differences for the case of A-AFM order. 

Fig.[4]shows the diagonal on-site energies e correspond- 
ing to Mn(e^), Mn{t2g), O(p^), and O(p^) MLWFs, as 
well as the on-site off-diagonal matrix elements q between 
the two Mn(e^) orbitals on the same site as function of 
the JT amplitude /Qq. It can be seen that the on-site 
energies for the various orbitals are essentially unaffected 
by the JT distortion. The corresponding changes, which 
lead to small differences between formerly symmetry- 
equivalent orbitals (see below for the case of the two Cg- 
like MLWFs), are of the order of a few tens of meV, and 
are thus negligible compared to the eg-t2g splitting or the 
d-p energy separation. 

A more distinct effect of the JT distortion can be 
seen in the on-site off-diagonal matrix element q be- 
tween the two majority spin Mn(e^) orbitals on the same 
site. This matrix element is zero for the undistorted 
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FIG. 5: Effect of increasing JT distortion on the hopping 
ampUtudes tt2g,p7r (a) and teg,po- (b) between Mn{d) and 0(p) 
nearest-neighbor orbitals (symbols). The different hopping 
ampUtudes are denoted via the corresponding Mn((i) orbital 
and the direction of the hopping ("in-plane" corresponds to 
x-y directions, "out-of-plane" to the ^-direction). The solid 
lines correspond to the hopping amplitudes for a perfect cubic 
perovskite structure with lattice constant a = ao =b Aa, where 
Aa = V2Q^ . Upper/lower branches correspond to short/long 
Mn-0 bonds. 



FIG. 4: On-site MLWF matrix elements of the d-p model as 
function of the JT distortion. (a)/(b): On-site energies for 
majority /minority spin, (c): On-site off-diagonal element q. 
(d): Splitting Ae between on-site energies of the two e^-like 
MLWFs. 



structure, but exhibits a linear increase to a value of 
about 0.25 eV for full JT distortion. This is consis- 
tent with the usual crystal-field picture of Cg orbitals in 
an octahedral environment (see also Eq. ([l])):^^ within a 
{|3z^— r^), \x'^ — y'^)} orbital basis, the Q^-type JT distor- 
tion gives rise to a nonzero off-diagonal matrix element 
that increases linearly with the JT distortion, while the 
diagonal elements of the Hamiltonian remain constant. 

However, the effect of the JT distortion observed 
for the minority spin q is significantly weaker than 
for the majority spin case. In fact, it is of simi- 
lar magnitude as the small splitting between the cor- 
responding diagonal matrix elements (on-site energies) 
Ae = £[Mn{3z'^ - r^)] - e[Mn{x'^ - y'^)]. We thus define 
a total JT-induced orbital splitting within the "e^" or- 
bital subspace as the difference in eigenvalues of the cor- 
responding 2x2 on-site Hamiltonian matrix, which is 
given hy S = \/ Ae'^ + Aq^. From the results presented 
in Fig. [4] for full JT distortion, i.e. Q^/Qg = 1, we ob- 
tain 5^ = 0.50 eV for majority spin and = 0.07 eV for 
minority spin. 

To understand this large difference between majority 
and minority spin, we note that in a partially covalent 
system like LaMnOa, the JT-induced orbital splitting 
is caused by a superposition of two effects: i) the true 
(electro-static) crystal- field effect, and ii) a "ligand- field" 
effect due to changes in hybridization with the surround- 
ing ligand orbitals. While the pure electro-static crystal- 
field effect is identical for both spin projections, the hy- 



bridization with the surrounding ligand orbitals is differ- 
ent for majority and minority spin orbitals, due to their 
large energy difference of about 3-4 eV. 

In order to verify whether the large difference be- 
tween majority and minority spin JT-splitting is indeed 
due to a strong ligand- field contribution, we have con- 
structed an alternative set of 20 MLWFs per formula unit 
and spin channel, where we explicitly included also the 
bands corresponding to 0(5) and semi-core La(5p) states. 
These states are energetically lower than the Mn{d) and 
0{p) bands shown in Fig. |2j and are centered around 
E ^ -3 eV (La(5p)) and E ^ -5 eV (0(s)). The in- 
clusion of these states in the MLWF basis reduces the 
"tails" in the e^-like MLWFs located at surrounding O 
and La atoms, and should therefore decrease the ligand- 
field contribution to the local JT splitting. Indeed, we 
obtain S^/^ = 0.36/0.12 eV in this case, i.e. the spin- 
dependent ligand-field contribution to the JT splitting 
is indeed significantly reduced compared to the d-p ML- 
WFs. The remaining difference can be ascribed to fur- 
ther contributions of orbitals that are not included in the 
MLWF basis and to the remaining small 0(5), 0(p), and 
La(p) contributions on the surrounding sites, which are 
required to ensure orthogonality between the MLWFs. 

We also note that the JT splitting corresponding to the 
rather localized Cg MLWFs of the d-p TB model (0.50 eV 
for majority spin) is a factor of two smaller than the 
corresponding splitting for the more spatially extended 
MLWFs of the effective Cg model which was reported in 
Ref. m (0.97 eV for majority spin). This demonstrates 
the much stronger ligand-field effect in the less localized 
effective Cg MLWFs. 

Next, we analyze the effect of the JT distortion on the 
nearest neighbor hopping between Mn{d) and 0{p) or- 
bitals. Fig. [5] shows the corresponding changes in the 



majority spin hopping amplitudes (the results for minor- 
ity spin do not exhibit any qualitative differences). 

We first discuss hopping between Mn{t2g) and 0(^7^) 
orbitals (Fig. [S^i). It can be seen that the hopping am- 
plitude tt2g,p7T splits essentially according to the different 
Mn-0 bond lengths in the JT distorted structure. The 
hopping amplitudes along the short (long) Mn-0 bonds 
within the x-y plane increase (decrease) with increasing 
JT distortion, while the hopping along the z-direction, 
i.e. corresponding to constant Mn-0 bond distance, is 
only weakly affected. 

In order to assess to what extent the d-p hoppings in 
the JT distorted structure are simply determined by the 
corresponding Mn-0 distances, we also calculate ML- 
WFs for undistorted cubic LaMnOs with different lat- 
tice constants a. The Mn-0 distances in the cubic struc- 
tures with a = ao ± Aa are identical to the long/short 
Mn-0 bond length in the JT distorted structure with 
= Aa/^/2. The corresponding results for the nearest 
neighbor hopping amplitudes between Mn{d) and 0{p) 
orbitals are also shown in Fig. |5j 

It is apparent that the tt2g,p7r-type hopping amplitudes 
in the JT distorted structures are nearly identical to the 
corresponding hopping amplitude in the undistorted cu- 
bic structure with the same Mn-0 distance. The small 
deviations between these two cases as well as the weak 
effect of the JT distortion on the tt2g,p7r hopping along z 
are due to changes in the orbital character of the MLWFs 
with increasing JT distortion. Thus, the magnitudes of 
the various hopping amplitudes are indeed determined 
mostly by the corresponding Mn-0 bond lengths. 

The case of the teg^pa hopping is only slightly dif- 
ferent. The hopping (both in-plane and along z) be- 
tween the Mn(32;^ — r^)-type orbital and the surrounding 
0{pa) orbitals compares well with the hopping ampli- 
tude in the undistorted structure with the same Mn-0 
distance (even though the agreement for the short Mn- 
O distance is not as good as for tt2g,p7r)- On the other 
hand, the JT-induced splitting of the in-plane hopping 
between the Mn(x^ — ?/^)-type orbital and the surround- 
ing 0{p) orbitals is significantly weaker than the corre- 
sponding bond-length dependence in the cubic structure. 
This indicates a strong change of the — y^)-type or- 
bital with increasing JT distortion. This change is due to 
different admixture of other orbitals that are not explic- 
itly included in the MLWF construction, i.e. 0{s) and 
La(p), which lead to a reduction/expansion of the lobes 
directed along the shorter/longer Mn-0 bonds with in- 
creasing JT distortion. This partially compensates the 
effect of changing Mn-0 distance and leads to the ob- 
served JT dependence of the hopping. 



C. Effect of U within the d-p model 

In order to investigate the influence of the Hubbard 
U parameter on the TB parameterization obtained from 
the MLWFs, we now perform GGA+/7 calculations for 
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FIG. 6: (Color online) On-site energies corresponding to the 
d-p model as a function of Hubbard U for majority spin (a) 
and minority spin (b). (c): U dependence of the on-site off- 
diagonal elements q. The small open (red) symbols repre- 
sent the GGA+[/ potential shifts added to the corresponding 
MLWF on-site matrix elements for U — 0. 



the fully JT distorted structure and FM order using dif- 
ferent values for U. Figure [6] shows the resulting U de- 
pendence of the on-site MLWF matrix elements. It can 
be seen that the diagonal elements (on-site energies e) 
corresponding to 0(p)-type MLWFs are virtually unaf- 
fected by the value of [/, while the corresponding matrix 
elements for Mn{t2g) and Mn(e^) orbitals exhibit a more 
or less linear dependence on U. 

This linear dependence is a direct consequence of the 
U dependent potential shift applied to the TM d states 
within the DFT approach: 

AC^. = t/(^-n^^O . (6) 

Here, n^^, is the occupation matrix element between 
atomic orbitals m and for spin projection a, which 
is calculated from the projection of the occupied Bloch 
states on fixed atomic orbitals. 

In order to quantify, whether the U dependence of the 
on-site MLWF matrix elements does indeed correspond 
to this potential shift, we explicitly evaluate Eq. (|6| for 
each [/, using the atomic occupation matrix elements ob- 
tained from the corresponding GGA+/7 calculation, and 
add the so-obtained shifts to the on-site MLWF matrix 
elements for U = eV. The resulting data is shown as 
small open (red) symbols in Fig. [6] and is nearly identical 
to the MLWF matrix elements calculated for the corre- 
sponding values of U. Thus, the GGA+/7 potential shifts 
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TABLE I: Average occupations corresponding to 

Mn(t2^) and Mn{eg) orbitals calculated for [/ = in the fully 
JT distorted structure. Rows denoted "atomic" correspond 
to the fixed atomic orbitals used to evaluate the GGA+?7 
functional; rows denoted "MLWF" contain the occupation of 
the MLWFs calculated according to Eq. ([7|); rows denoted 
"formal" correspond to the ionic limit based on a high-spin 
configuration of the Mn cation. The last column contains 
the off-diagonal occupation matrix element n^^/ between the 
two different Cg orbitals. The first/last three rows correspond 
to majority /minority spin. 
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are directly reflected in the on-site matrix elements of the 
MLWFs. 

The reason for this good correspondence between the 
GGA+/7 potential shifts and the U dependence of the on- 
site MLWF matrix elements is the fact that the MLWFs 
of the extended d-p model are rather similar to the atomic 
orbitals used as projector functions within the GGA+/7 
approach. To further demonstrate this similarity, we also 
compare the occupation matrix elements used to evaluate 
the GGA+/7 potential shifts with the occupation of the 
corresponding MLWFs 



MLWF 



BZ 7 



(7) 

where is the Fermi energy. The corresponding values 
for U = eV are listed in Table |Tj It is apparent that 
the occupations of the MLWFs are very similar to the 
occupations of the atomic orbitals used as GGA+/7 pro- 
jector functions. Furthermore, we note that due to the 
strong hybridization between the Mn{d) and 0{p) states, 
the occupation of the atomic-like d orbitals are quite dif- 
ferent from a naive expectation based on the formal ionic 
configuration of the Mn^+ cation. 

The U dependence of the on-site off-diagonal matrix 
element q can be understood in the same way and is 
related to the off-diagonal occupation matrix element 
^mm'' where m=\3z'^ — r^) and m' = \x'^ — y'^). The 
nonlinearities that can be observed for large U are due 
to changes in the corresponding n^^,. Note that since 
q is directly related to the JT splitting, increasing U ef- 
fectively amounts to increasing the strength of the JT 
coupling. 

The U dependence of the hopping amplitudes between 
0{p) and Mn((i)-like MLWFs for majority spin and FM 
order is shown in Fig. [7[ All depicted hopping ampli- 
tudes exhibit a similar small increase in magnitude. The 
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FIG. 7: Magnitude of hopping amplitudes between Mn(d) and 
0(p) nearest-neighbor orbitals as function of the Hubbard U. 
(a): t2g —p-K hoppings, (b): Cg —pa hoppings. Notation is the 
same as in Fig. |5] 



effect of U on the hopping amplitudes for minority spin 
(not shown) is even weaker than for the majority spin 
orbitals. These small changes in the hopping amplitudes 
for both majority and minority spin can be attributed 
to changes in orbital character of the corresponding ML- 
WFs. These changes in orbital character arise from a 
different admixture (hybridization) of orbitals centered 
at the surrounding ions, which results from the U de- 
pendent shifts of the on-site energies of the Mn{d) states 
relative to all other orbitals. 



D. Effect of U within the effective Cg model 

We now contrast the U dependence of the TB parame- 
ters for the extended d-p model presented in the previous 
section with the case of an effective Cg model with only 
two e^-like Wannier orbitals on each Mn site. This is 
probably the most common model used to describe the 
physics of manganites (see e.g. Ref. [3|, and a detailed 
analysis of the effect of various distortions and different 
magnetic arrangements on the corresponding MLWF ma- 
trix elements for U = eV has been presented in Ref. |9l 
Here we focus on changes in the MLWF parameteriza- 
tion due to the Hubbard U correction for the purely JT 
distorted structure and FM order. The corresponding 
Wannier functions are constructed from the Kohn-Sham 
states located within an energy window of 12.0-17.0 eV 
and 15.9-20.0 eV for majority and minority spin, respec- 
tively. For further details see Ref. 9, 

Fig. [S] shows the U dependence of the on-site MLWF 
Hamiltonian matrix elements for the effective Cg model. 
It can be seen that in particular the on-site energies ex- 
hibit a very similar trend as the on-site energies of the 
e^-like MLWFs in the d-p model, but that the U de- 
pendence is weaker than in the latter case. The slope 
de/dU for both majority and minority spin is only about 
66 % compared to the d-p model. This indicates that 
the U dependence of the on site energies for the effective 
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(a) (b) 




U (eV) U (eV) 



FIG. 8: On-site energies e (a) and on-site off-diagonal matrix 
element q (b) of the effective Cg MLWF parameterization as 
function of the Hubbard U . Majority spin and minority spin 
is denoted by up and down triangles, respectively. 



Cg MLWFs is still determined by the GGA+[/ potential 
shifts, Eq. (|6|, but is renormalized by the extent of over- 
lap between the extended MLWF and the corresponding 
eg atomic orbital. In other words, the U dependent en- 
ergy shift experienced by the MLWF is determined by 
the projection of the MLWF on atomic orbit als and the 
occupation of these atomic orbitals. 

In particular, this shows that the effect of U on the 
on-site energies of the effective eg MLWFs are not de- 
termined by the occupation of the MLWFs themselves 
via a relation similar to Eq. ([6|. The occupation of the 
effective eg MLWFs are essentially identical to the "for- 
mal" occupations listed in Table |Tj corresponding to the 
formal 3+ charge state of the Mn cation in LaMnOa (0.5 
for majority spin and 0.0 for minority spin). Therefore, 
the occupation-dependent potential shift following from 
a mean-field approximation to the Hubbard interaction 
similar to Eq. (|6| would be zero for the majority spin 
eg states. The U dependence of the corresponding on- 
site energies is thus notably different from a mean-field 
treatment of the electron-electron interaction in a two- 
orbital Hubbard-like model derived from the effective eg 
TB parameterization. 

We point out that in addition to the different de- 
pendence on orbital occupation, the magnitude of the 
screened Hubbard U acting on the extended e^-like Wan- 
nier orbitals will of course be significantly reduced com- 
pared to the magnitude of U corresponding to more lo- 
calized atomic-type orbitals. A comparison of the value 
of U for "atomic-like" and "effective" Wannier orbitals 
in an Fe-pnictide system based on calculations using con- 
strained random phase approximation has been presented 
recently in Ref. 35, 

Similar to the case of the on-site energies, the U de- 
pendence of the off-diagonal on-site matrix elements 
shown in Fig. (8)3, resembles the U dependence of the cor- 
responding atomic-like MLWFs shown in Fig. [Bj^. The 
decrease of q for majority spin observed ioi U > 6 eV 
can then be ascribed to changes in orbital character of 
the extended MLWFs, which lead to different projections 
between MLWFs and the corresponding atomic orbitals. 




U (eV) U (eV) 

FIG. 9: Effect of U on the nearest neighbor hopping ampli- 
tudes of the effective Cg MLWFs for majority spin (a) and 
minority spin (b) . "In-plane" and "out-of-plane" corresponds 
to hopping between similar orbitals within the x-y plane and 
along ^, respectively, while "in-plane mixed" denotes the hop- 
ping between two different Cg MLWFs at neighboring sites 
within the x-y plane. 



Finally, Fig. [9] shows the U dependence of the near- 
est neighbor hopping parameters for the effective eg ML- 
WFs. For majority spin, there is an overall increase in 
the magnitude of the hopping, whereas for minority spin 
there is an overall decrease. Compared to the d-p model 
discussed in the previous section, the effect of U on the 
hopping amplitudes of the effective eg model is strongly 
enhanced, with typical changes between 30 % and 60 % 
over the considered range of U values. For the in-plane 
hopping between |3z^ — r^) and — ?/^)-type orbitals 
the change can be even more than 100%. 

This strong effect of U on the hopping amplitudes of 
the effective eg MLWFs can be explained by the changes 
in the underlying d-p model. The effective eg MLWFs are 
essentially linear combinations of an atomic-like Mn(e^) 
orbital and O(pcr) orbitals on the surrounding anions. As 
demonstrated in the previous section, the on-site energies 
of the atomic d orbitals are strongly shifted relative to 
the oxygen p orbitals as a function of U . This shift af- 
fects the hybridization between the atomic orbitals and 
thus changes the admixture of O(pcr) orbitals in the effec- 
tive eg MLWF. For majority spin, increasing U decreases 
the energy difference between atomic 0{p) and Mn((i) 
orbitals and thus increases the hybridization, i.e. the ad- 
mixture of 0{p) in the effective eg MLWF. This leads to 
an increase of the hopping amplitude due to the larger 
overlap between effective eg orbitals at neighboring Mn 
sites. For minority spin the trend is opposite to this, 
resulting in the decrease in magnitude seen in Fig. ^p. 

Considering a simple d-p nearest neighbor TB model 
for the cubic perovskite structure, one can show that in 
the limit of large energy separation between the d and p 
orbitals, the "antibonding" bands that result from the d-p 
hybridization are formally equivalent to an effective "(i- 
only" TB model with direct hopping between d-orbitals 
on adjacent TM sites. The amplitude of this "effective" 
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d-d hopping is thereby given as: 



Obviously, the overah trends discussed in the previous 
paragraph are consistent with this simple picture, but 
the limit for which Eq. (|8| is valid is not fulfilled in the 
case of LaMnOa. This is apparent from the on-site ener- 
gies shown in Fig. [6^, where it can be seen that for values 
of U larger than 5-6 eV the Mn(e^) states become essen- 
tially degenerate with the 0{p) states. We also verified 
that Eq. ([8| is not fulfilled quantitatively by the hoppings 
calculated from both sets of MLWFs {d-p and effective 
eg). Thus, even though it is possible to obtain a reason- 
able MLWF parametrization of the effective Cg bands for 
all values of /7, it is not obvious whether a low-energy 
description of LaMnOs based only on effective Cg states 
is always physically reasonable. 

The partial "breakdown" of the effective Cg description 
of LaMnOs for values of U larger than 5-6 eV also leads to 
a strong dependence of some MLWF parameters on sub- 
tle differences in the underlying Kohn-Sham bandstruc- 
ture. For example the "crossing" of the mixed in-plane 
hopping amplitudes, i.e. corresponding to the in-plane 
hopping between |3z^ — r^) and — y^)-type orbitals, 
that can be seen around U = 5 eV in Fig. [9^, appears 
only if orthogonalized atomic orbital are used to evaluate 
the "-j-t/" correction to the GGA energy functional. If 
instead a projection on non-orthogonal atomic orbitals is 
used (option "atomic" instead of "ortho-atomic" within 
QuantumESPRESSO), the resulting MLWF matrix ele- 
ments do not exhibit this feature. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have discussed differences in the 
MLWF-derived TB parameterization of LaMnOs that 
result from different values of the Hubbard U used in 
the GGA+/7 calculation from which the MLWFs are ob- 
tained. Thereby, we have compared two different ways 
to represent the important bands around the Fermi en- 
ergy. First, a d-p TB model based on atomic-like Mn{d) 
and 0{p) MLWFs, and second, an effective Cg TB model 
involving only two e^-like Wannier orbitals per Mn. 

We have shown that the hopping amplitudes of the 
d-p model are only weakly affected by a variation of /7, 
and that the resulting changes in the on-site Hamiltonian 
matrix elements are consistent with a mean-field approx- 
imation to the electron-electron interaction for the cor- 
responding TB Hubbard Hamiltonian. As a result, the 
TB parameters for the d-p model of LaMnOs are fairly 
insensitive to variations of U and can therefore be deter- 
mined accurately without detailed knowledge about the 
precise value of U for the corresponding d orbitals. 

In contrast, the hopping amplitudes for the effective 
Cg TB parameterization depend strongly on the value 



of U used in the GGA-\-U calculation. This is due to 
pronounced changes in the amount of admixture of 0{p) 
ligand orbital character in the effective Cg MLWFs, which 
in turn is due to the large U dependent shifts of the on- 
site energies of the corresponding atomic orbitals. Fur- 
thermore, the U dependence within a mean-field Hub- 
bard Hamiltonian derived from the effective Cg MLWF 
parameters is distinctly different from the U dependence 
of the MLWF parameters themselves. The former is de- 
termined by the occupations of the effective Cg orbitals 
and the value of the screened Hubbard interaction in that 
basis, whereas the latter is determined by the occupa- 
tions and the value of the screened Hubbard interaction 
corresponding to atomic-like d states. A determination 
of suitable TB parameters for an effective Cg model of 
LaMnOs therefore requires an accurate knowledge of U 
in both basis sets. This is particularly important for the 
determination of the JT coupling strength, which is de- 
termined from the on-site splitting within the Cg orbital 
manifold. This splitting on the other hand is critically af- 
fected by the value of U in both the GGA+/7 calculation 
and the mean-field Hubbard Hamiltonian. 

In addition to the U dependence, we have also an- 
alyzed the effect of the staggered JT distortion ("Q^- 
type" ) on the TB parameterization of the d-p model. We 
have found that the JT distortion manifests itself both 
as a local ligand-/crystal-field splitting within the Cg or- 
bital manifold, as well as via a pronounced bond-length 
dependence of the nearest neighbor d-p hopping ampli- 
tudes. Our results demonstrate that even for the rather 
localized atomic-like d-p MLWFs, the ligand-field effect 
dominates over the purely electro-static crystal-field ef- 
fect. Furthermore, we have verified that the changes in 
the nearest neighbor hopping amplitudes due to the JT 
distortion are almost fully determined by the resulting 
change in Mn-0 bond length. 

A number of more general conclusions can be drawn 
from the results obtained in this work. It is reason- 
able to assume that the sensitivity /insensitivity of cer- 
tain MLWF parameters from the specific value of U used 
in the GGA+/7 calculation can be generalized to a gen- 
eral sensitivity /insensitivity from the choice of exchange- 
correlation functional used in the electronic structure cal- 
culation. In the present study, we made use of the sim- 
ilarity between the "+/7" correction to the GGA func- 
tional and the mean-field approximation to the electron- 
electron interaction in the corresponding model Hamil- 
tonian (assuming of course that the latter is described 
via a local Hubbard interaction). Using this similarity, 
the observed trends in the MLWF TB parameters, and 
the difference between the two different TB models can 
be explained. The same similarity between DFT and 
mean-field Hubbard model is not present if the electronic 
structure is calculated using other "beyond LDA/GGA" 
methods such as for example self-interaction corrected 
methods, hybrid functionals,!^ or the GW approach.!^ 
In these cases the electronic band-structure should also 
be compared with the mean-field approximation of an 
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FIG. 10: (Color online) Schematic depiction of all closest 
neighbor Mn-0 and 0-0 hopping amplitudes within the ex- 
tended d-p model of LaMnOs- 



appropriate TB Hamiltonian. To extract the "bare" or 
"non-interacting" on-site energies and off-diagonal on- 
site matrix elements, one has to then subtract a U de- 
pendent shift similar to Eq. [6j The corresponding value 
of U has to be calculated by other means, for example 
via the recently proposed approach based on constrained 
random phase approximation.^^ 

While the so-obtained TB models will always give 
an essentially perfect representation of the Kohn-Sham 
bands for the given reference, it is still important to 
test whether the corresponding parameterization is also 
transferable to slightly different configurations, i.e. with 
different structural distortions and/or magnetic/orbital 
order. Of course, as shown through the comparison be- 
tween the effective Cg and the more elaborate d-p model 
for LaMnOs, transferability can always be improved by 
including more orbit als in the TB basis set. 

Finally, we note that in the present work we have fo- 
cused on how to obtain TB parameters using the DFT 
Kohn-Sham band-structure as (mean-field) reference. An 
entirely different, albeit at least equally important, ques- 
tion is whether these bands, calculated using a suitable 
exchange-correlation functional for a particular system, 
are indeed a good representation of the real material. 
However, this question is beyond the scope of the present 
study and will most likely be the topic of future research 
for a number of years to come. 
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Appendix A: Further details on the d-p TB 
parameterization for cubic LaMnOs 

In this appendix we give a more quantitative account 
of the d-p TB parametrization obtained from the MLWFs 



TABLE II: On-site energies and nearest neighbor Mn-0 and 
0-0 hopping amplitudes obtained for the d-p MLWFs in the 
cubic structure with FM order. In addition, the hopping be- 
tween t2g orbitals along the shortest Mn-Mn distance is also 
listed. All values are given in eV. 
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FIG. 11: (Color online) Band dispersion calculated from the 
MLWFs (thick line) and from the truncated TB model de- 
scribed in the main text (dots and line); (a): majority spin, 
(b): minority spin. 



for simple cubic LaMnOs and FM order. 

There are 2 different hopping parameters between 
nearest neighbor Mn{d) and 0{p) states {teg, pa ^ ^t2^,p7r) 
and 4 independent hopping parameters between closest 



neighbor 0{p) states (t 



t 



pa, pa-) ^pa^pTT") ^pTv^pTV) ^P7v,p7r 



ti 



that 



are allowed by the symmetry of the system. All of these 
hopping parameters are illustrated in Fig. [lO) The corre- 
sponding values together with the various on-site energies 
are summarized in Table Ull It can be seen that the mi- 
nority spin hopping amplitudes corresponding to nearest 
neighbor d-p hopping are slightly larger than for major- 
ity spin, whereas the opposite is the case for the closest 
neighbor p-p hoppings (except for t^j^^pT:). 

Within the two-center approximation to the linear 
combination of atomic orbitals (LCAO) method all clos- 
est neighbor 0{p)-0{p) hoppings can be expressed as lin- 
ear combination of purely a- and 7r-type hopping param- 
eters, t{ppa) and t{pp7r).^^ Thereby, t{pp7r) corresponds 
directly to iIt^^ptt in our notation. Calculating an esti- 
mate for t{pp(j) from each of the remaining MLWF hop- 



ping parameters as (t 



pTTjPTT ' 



it. 
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-2ti 



r). or 
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{—tp'K,p'K + '^tpcr^pTv), we obtain values of 0.873/0.725 eV, 
0.660/0.516 eV,'or 0.627/0.622 eV, for majority /minority 
spin, respectively. The spread in these values shows to 
what extent the approximation of rigid atomic orbitals is 
not fulfilled within the d-p set of MLWFs. 

In a previous study, LCAO parameters for manganites 
were deduced based on a cluster-model analysis of pho- 
toemission spectra,^^ leading to the following values for 
the hopping integrals: t{pda) = 1.8 eV, t{pd7r) = 0.9 eV, 
t{ppa) = 0.60 eV, t{pp7r) = —0.15 eV. These values are 
in reasonable agreement with the values obtained from 
the MLWFs listed in Table III 



In Figure 10 we show a comparison of the dispersion 
between the MLWF bands and a 14-band TB model cal- 



culated using the parameters given in Table |Tl| i.e. by 
neglecting all other further neighbor hoppings. The over- 
all band dispersion is rather well reproduced in the trun- 
cated TB model. However, due to a rather slow decay of 
hopping amplitudes with further neighbor distance, cer- 
tain features are not reproduced. In addition, neglecting 
further-neighbor hoppings also leads to a slight under- 
estimation of the total bandwidth in the truncated TB 
model. We find that there are overall 14 different fur- 
ther neighbor hoppings with magnitudes in the range of 
(0.020-0.113) eV. Out of these only the direct hopping 
between neighboring t2g orbitals (last line in Table M 
leads to a significant improvement of the TB bands and 
is therefore included in Table [ill 
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would be different, our general conclusion would be un- 
affected. For simplicity of presentation we therefore only 
discuss the case J = 0. 
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